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ABSTRACT 

We report a large spin-up glitch in PSR J1846— 0258 which coincided with the onset of 
magnet ar- like behavior on 2006 May 31. We show that the pulsar experienced an unusually 
large glitch recovery, with a recovery fraction of(3 = 8.7±2.5, resulting in a net decrease of the 
pulse frequency. Such a glitch recovery has never before been observed in a rotation-powered pul- 
sar, however, similar but smaller glitch over-recovery has been recently reported in the magnetar 
AXP 4U 0142-f 61 and may have occurred in the SGR 1900-1-14. We also report a large increase 
in the timing noise of the source. We discuss the implications of the unusual timing behavior in 
PSR J1846— 0258 on its status as the first identified magnetically active rotation-powered pulsar. 

Subject headings: 



pulsars: general — pulsars: individual (PSR J1846— 0258) — X-rays: stars 



Introduction 



PSR J1846-0258 is a young ('-800yr), 326 ms 
pulsar, discovered in 2000 with the Rossi X-ra y 
Timing Explorer {RXTE] iGotthelf et all |2000[) . 
No radio puls ations have been detected despite 
deep searches ( Kaspi et al.lll996t Archibald et al. 
20081 ). PSR J1846-0258 has a large inferred 
magnetic field of i? ~ 5 x 10^^ G, above the 
quantum critica l limit, and a braking index of 
n = 2.65 ±0.01 (jLivingstone et al.ll2006l hereafter 
LKGK06). The pulsar was therefore believed to 
be simply a high-magnetic field rotation-powered 
pulsar (RPP), with radio pulsations that do not 
cross our line of sight. However, in 2006 May, the 
source experienced a series of X-ray bursts and a 
sudden increase in X-ra y flux, which strong ly sug- 
gest magnetic activity ( Gavriil et al.l 120081 ). The 



source also showed spectral va riations and changes 
in the surrounding nebu la (iGavriil et al. 2008t 
Kumar fc Safi-Harbll2008HNg et al.ll2008l) . 



A magnetar is a 
ne tic energy powers 



neutron star whose mag- 
the bulk of its emission 
while a RPP pro- 
of rotational ki- 



(jThompson fc DuncanI 119951) . 

duces radiation vi a the loss 

netic energy (e.g., Manchester fc Tavloil Il977() . 
PSR J1846— 0258 is a unique transition object 
between these two source classes. The X-ray lu- 
minosity of PSR J1846— 0258 can be entirely ac- 
counted for by the spin-down power of the pulsar, 
however, the observed X-ray bursts and flux flare 
are phenomena only seen thus far from magnetars 



(|Gavriil et al.ll2008D . 
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A neutron star glitch is defined as a sudden, 
usually unresolved, increase in spin frequency, v. 
Glitches are often accompanied by an increase in 
the frequency derivative, iy, and are sometimes fol- 
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lowed by an exponential decay on timescale Td, 
where some of the initial jump in v is recovered. 
In general, a glitch at time tg can be modeled as: 

-(t-ts)/Td 



(1) 

where vo(t) is the frequency of the pulsar prior to 
the glitch and Av is the initial frequency jump, 
which can be decomposed into the part of the 
glitch that is permanent, Avp, and that which 
decays, Az/^- The recovery fraction is defined 
as Q = Avdl Av. The variety of observed pul- 
sar glitches is interesting. The fractional mag- 
nitude of the observed change in v ranges over 
6 ord ers of magnitude from 10~^^ < Avjv < 
10"^ (iLvne et alJ[2000Hjanssen fc Stappers|[200l 
Hobbs et al.l 2002 ). Some glitches are character- 
ized only by a ly incr ease (e.g., PSR B1758— 23; 



IShemar fc L^ll996h . while some glitches (typi- 
cally in young pulsars) are dominated by an in- 
crea se in the magnitud e of v (e.g., the Crab pul- 
sar; IWong et"aL I l200lh . S ome glitches de cay on 
timescales of '^days (e.g., iFlanagaii 199Cll ). while 
others display very long decay tim e scales (~ hun- 
dreds of days; IWang et aL I I2OOOI) . The fraction of 
the glitch that decays is also highly variable. In 
older pulsars the amo unt of glitch recovery i s typi- 
cally small (Q « 1, IShemar fc Lvnel[l99fih while 
in the Crab pulsar the recovery can b e nearly com- 
plete (e.g., Q ^ 0.96, iLohsenlllQSlI ). However, in 



many glitches, no recovery is detected. 

Glitches are now known to be ubiquitous 
in m a.gnetars as well a s rotation-powered pul- 
sars (iKaspi e t al.' '2000"; 'Kasoi fc Gavriil' l2003t 
DairOsso et aL 2003: Dib et al. 2008: Gavriil et al. 



20091 ). What remains to be seen is whether 
glitches have the same physical origin in both 
types of objects, or if the super-critical magnetic 
fields of magnetars are responsible for different 
glitch origins and evolutions. While some mag- 
netar glitches are indistinguishable from those 
observed in RPPs, others occur contemporane- 
ously with radiative changes such as bursts, flux 
enhancements, and spectral or pulse profile vari- 
ations, such as in th e Anomalous X-ra y Pulsar 
(AXP) IE 2259-^586 (jKaspi et al.ll2003[) . No ra- 
diative chang es have been obser ved with RPP 
glitches (e.g., iHelfand et alll200l[ ). RPP glitches 
are believed to arise from a sudden unpinning 
of vortices in t he superfluid interior crust of 
the pulsar (e.g., lAlpar et al.l 119931 ). Magnetar 



glitches, on the other hand, may be triggered 
by strong internal magnetic fields as the crust 
i s deformed, either plasticall y or cracked violently 
(jThompson fc Duncanlll996l ). 

In this paper, we discuss the timing behavior of 
PSR J1846— 0258 prior to, during, and following 
the period of magnetic activity observed in 2006. 
We show that a large glitch occurred contempora- 
neous with the X-ray bursts and onset of the flux 
flare, accompanied by an unusual increase in the 
timing noise of the pulsar. We show that the glitch 
recovery is very unusual for a RPP but is remi- 
niscent of timing behavior observed from magne- 
tars and is further evidence of magnetic activity 
in PSR J1846-0258. 

2. RXTE Observations and Analysis 

Observations of PSR J1846-0258 were made 
using the Propor t ional Counter Arr ay (PGA; 
Jahoda et al. 19961 Jahoda et al.ll2006l) on board 
RXTE. The PGA consists of an array of five col- 
limated xenon/methane multi-anode proportional 
counter units (PCUs) operating in the 2 - 60keV 
range, with a total effective area of approximately 
6500 cm^ and a field of view of - 1° FWHM. 

Our entire RXTE data set spans 9.7 yr from 
1999 April 18 through 2008 December 10 (MJD 
51286 - 54810). Data from 2000 January 31 - 
2005 July 27 (MJD 51574 - 53578) were reduced 
and analyzed previously and details can be found 
in LKGK06. Eleven observations taken in 1999 
April 18-21 are of limited use for the current anal- 
ysis since they cannot be unambiguously phase 
connected to the rest of the data. Analysis of 
data spanning 2005 July 27 - 2008 December 10 
(MJD 53578 - 54810) is presented here. Data were 
collected in "GoodXenon" mode, which records 
the arrival time (with l-/is resolution) and en- 
ergy (256 channel resolution) of every unrejected 
event. Typically, two to three PCUs were op- 
erational during an observation. We used the 
first Xenon layer of each operational PGU and ex- 
tracted events in channels 4-48 (approximately 
2 - 20kev), as this maximizes the signal-to-noise 
ratio for this source. 

Observations were downloaded from the HEASARC 
archiv£0 and data from each active PGU were 



^http: / /heasarc. gsfc.nasa.gov/docs/archive. html 
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merged and binned at (1/1024) s resolution. Pho- 
ton arrival times were converted to barycen- 
tric dynamical time (TDB) at the solar system 
barycenter using the J2000 source position RA = 
18M6'°24!94 lb 0!01 Decl = -02°58'30.1" ± 0.2" 
(jHelfand et all l2003l) and the JPL DE200 solar 



system ephemeris. 

The known ephemeris from LKGK06 was used 
to fold each time series with 16 phase bins. Re- 
sulting profiles were cross-correlated in the Fourier 
domain with a high signal-to-noise ratio template 
created by adding phase-aligned profiles from all 
observations. The cross-correlation process as- 
sumes that the pulse profile is stable; indeed, 
we found no evidence for variability that could 
bias TOA measu rement, as shown i n Fig ure [T] 
and confirmed bv lKuiper fc HermsenI (120091 ). The 
figure shows average profiles from before, after, 
and throughout the outburst and unusual tim- 
ing behavior. For each individual observation, the 
cross-correlation yielded the time of arrival (TOA) 
of phase-zero of the average pulse profile at the 
fold epoch. The TOAs were fitted to a timing 
model (see Section [S]) using the pulsar timing soft- 
ware package TEMPCd. After phase-connecting 
the data, we merged observations occurring on a 
single day and used the ephemeris to re-fold the 
data in order to obtain more precise TOAs. This 
process produced 199 TOAs with a typical uncer- 
tainty ^9 ms (■~2.7% of the pulse period). Further 
details of the observation and analysis are given in 
LKGK06. 

3. Timing Analysis and Results 

Phase-coherent timing is a powerful method 
for obtaining accurate pulsar parameters, but can 
only be used when t iming noise and glitche s are 
relatively small (e.g., Livingstone et al.l 20051) . To 
obtain a phase-coherent timing solution, each 
turn of the pulsar is accounted for by fitting 
TOAs with a Taylor expansion of the pulse phase 
()Lvne fc Smithll2005l) . 

Phase-coherent timing for PSR J1846-0258 
spanning 2000 January 1 - 2005 July 27 (MJD 51574- 
53578) is discussed in LKGK06. Phase coherence 
was maintained for the next 308 days without 
incident. Phase coherence was lost with the ob- 
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Pulse Phase (2 cycles) 

Fig. 1. — Five average 2-20 keV pulse profiles 
for PSR J1846-0258. Each profile is based on a 
single phase-coherent subset of data, except for 
the profile made from un-connected periodogram 
measurements of v after the burst. For each pro- 
file, the center date of observations is listed in the 
top right hand corner. The top panel shows the 
added profile for a phase-coherent timing solution 
spanning MJD 53228-53879, with i^, and i> fit- 
ted. The profile includes 29 observations with a 
total exposure time of 223 ks. The second panel 
shows the summed profile for 9 observations af- 
ter the X-ray bursts and flux flare spanning MJD 
53949-54005 with an exposure time of 44.3 ks. The 
third panel shows the profile for a phase-coherent 
timing solution of 8 observations with v and fit- 
ted spanning MJD 54215-54265, near the unusual 
variations in i/. The total exposure time is 38.9 ks. 
The fourth panel shows the average pulse profile 
near the end of the glitch recovery, using a coher- 
ent timing solution of 9 observations with and z> 
fitted, spanning MJD 54363 - 54397. The exposure 
time is 40.2 ks. The bottom profile shows the aver- 
age pulse profile after the glitch had recovered us- 
ing a coherent timing solution of 12 observations, 
spanning MJD 54726-54805. The exposure time is 
59.0 ks. 



^http: / /www. atnf.csiro.au/research/pulsar/tempo 
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servation occurring on 2006 May 31 (MJD 53886) 
which co ntained 4 X-ray bu rsts and a pulse flux 



increase ( Gavriil et al. 2008[ ). For the following 32 



observations spanning 192 days, no unambiguous 
phase coherent timing solution was possible. In- 
stead, we performed periodograms to determine 
the pulse frequency. Uncertainties were deter- 
mined from a Monte Carlo simulation, where noise 
was added to simulated sinusoidal pulses and the 
frequency for each trial was determined in the 
same way as for the real data (see LKGK06 for 
further details). Phase coherence was once again 
obtained starting 2007 January 26 (MJD 54126) 
with closely spaced bootstrapping observations 
after the source reappeared from behind the Sun 
after 48 days. A single phase coherent timing 
solution was obtained spanning 2007 January 26 
- 2008 December 10 (MJD 54126-54810). This 
timing solution is severely contaminated by long- 
term glitch recovery and timing noise, so the fitted 
parameters for the global solution are of limited 
value. 

To analyze the long-term rotational behavior of 
the pulsar, we created short phase-coherent tim- 
ing solutions from 2000 until the onset of bursts, 
and from 2007 and 2008. Each timing solution in- 
cluded only V and and included as much data 
as possible while requiring the reduced value of 
the fit to be ~1. This resulted in 11 measurements 
of V and v pre-glitch and 11 measurements post- 
glitch. In order to better utilize the available data 
in the post-glitch period, we created short over- 
lapping timing solutions, each of which uses ap- 
proximately half the data from two of the above 
described short data subsets, and has the same 
fitted parameters and requirements. This can 
be useful because the short coherent timing solu- 
tions result in parameter fits that are dominated 
by the end points, which can be problematic when 
V is varying rapidly from glitch recovery as in this 
case, or when timing noise is a significant effect. 
This produces an additional 9 post-glitch measure- 
ments of V and v. Coherent frequency measure- 
ments (crosses), overlapping frequency measure- 
ments (filled circles) and frequency measurements 
obtained via periodograms (open circles), are plot- 
ted in the top panel of Figured! with the pre-burst 
J/, !/, and V removed. The middle panel of Figure [2] 
shows measurements of v from the short coherent 
timing solutions as crosses, with the overlapping 



ii measurements in filled circles. In addition, open 
circles show three measurements of v that are cal- 
culated from weighted least-squares fits to nine pe- 
riodogram measurements of v. 

A frequency increase is apparent in the fre- 
quency residual plot (top panel and inset. Fig- 
ure [2|), i ndicating that a glitch o ccurred (also 
noted by Kuiper fc Hermsen 20091 ). Two mea- 
surements of V are larger than the pre-glitch pre- 
dicted value. We calculated the average I\i> = 
(5.5 ± 2.3) X 10-6 at MJD 53890 (4- 11 days after 
the glitch occurred). The initial value of Az/ at 
the time of the glitch was presumably larger than 
this value, as indicated by the exponential fit to 
the data discussed below. 

Wc fitted the measured frequencies of PSR J1846— 0258 
(fitting only 3 pre-glitch v values) with an expo- 
nential recovery glitch model (Eq. 1), the results 
of which are shown in Figure [3l The top panel 
of the Figure shows frequency measurements with 
the pre-glitch ephemeris {v, v and v) removed (as 
in Figure [2]), while the bottom panel shows resid- 
uals from the glitch fit. The uncertainties from 
the periodogram measurements of v are '^2 — 3 
orders of magnitude larger than those from short 
coherent fits to the data and thus contribute min- 
imally to the overall value. Thus the glitch 
residuals (bottom panel. Figure [3]) are shown to 
highlight the deviation from the fit of the coherent 
frequencies, which dominate the x^- 

Significant deviations from the fit can be seen 
during the period of glitch recovery, giving rise to 
a large reduced x^ value of ~ 267 for 45 degrees of 
freedom for the best fit (bottom panel. Figure [3l). 

Given the poor fit, the formal uncertainties un- 
derestimate the true values. To determine more 
reasonable uncertainties on the glitch parameters, 
we multiplied the uncertainties on the coherent v 
measurements by a factor until the reduced x^ of 
the fit ^ 1. We applied a multiplicative factor 
only to the phase coherent v measurements be- 
cause these uncertainties were ~ 2 — 3 orders of 
magnitude smaller than those from periodogram 
V measurements, determined in a different way, 
and the contribution to x^ was mainly from the 
coherent v values. Quoted uncertainties on the 
glitch parameters are from Ax^ = 1 contours 
from the fit with multiplied uncertainties. The 
fitted value of the initial fractional frequency in- 
crease is /S.v/v = 4.0(1.3) X 10^6, very large for 
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Fig. 2. — Timing and pulsed flux evolution of 
PSR J1846-0258 over 9 years. Top panel: Fre- 
quency measurements with pre-glitch v, and v 
removed from all points for clarity. Black crosses 
are produced from short phase-coherent measure- 
ments of V and while filled red circles are pro- 
duced similarly but from overlapping segments of 
data. Uncertainties are smaller than the points 
and are excluded for clarity. Open blue circles are 
periodogram measurements of v. Middle panel: 
Frequency derivative {u) measurements. Crosses 
are produced from short phase-coherent measure- 
ments of V and i>, while filled red circles are from 
overlapping segments of data, but are otherwise 
produced in the same manner. See Section [3] for 
details. Open blue circles are from weighted least- 
squares fits of periodogram v measurements. Bot- 
tom panel: Pulsed intensity measurements in the 
2 - 20keV energy band, as described in Section 
m The only significant flux increase is coincident 
with 4 X-ray bursts observed on MJD 53886. All 
panels: The dashed line represents the epoch of 
a small Crab-like glitch near MJD 52210, while 
the dotted line represents the epoch of a small 
candidate glitch between MJDs 52837 and 52915 
(LKGK06). The two solid l ines show t he epochs 
where bursts were detected ( Gavriil et al..i2008il . 
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Fig. 3. — Top panel: Frequency measurements of 
PSR J1846-0258 (points) and the fitted exponen- 
tial recovery glitch model (solid line). Pre-glitch 
V, and v have been removed from all frequency 
measurements. Bottom panel: Residuals from the 
fit with the formal, un-multiplied uncertainties. 
The best fit has a. ^ oi ^ 267 for 45 degrees 
of freedom. Significant variation from the fitted 
exponential is clear, however, the exponential re- 
covery dominates the change in by ~ 2 orders 
of magnitude over the remaining systematic varia- 
tions in V. Given the large , we increased the un- 
certainties on the phase-coherent v measurements 
until u ^ 1; Slid report uncertainties on the fit- 
ted model from the Ax^ = 1 contours. All fitted 
parameters are given in Table [TJ 
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such a young pulsar, where a typical value is 
/S.v/v ~ 10~®. More remarkable yet, however, is 
the amount by which the frequency recovers. We 
find Q — 8.7(2.5), corresponding to a net decrease 
in frequency of Az/p = 9.52(9) x 10"^ Hz. All fitted 
glitch parameters are given in Table [TJ 

While the deviation from the exponential fit is 
very significant for several months (bottom panel. 
Figure [3]) the overall evolution after the glitch is 
dominated by the exponential recovery: the de- 
viation from the fit is ^2 orders of magnitude 
smaller than the overall post-glitch decrease in v. 
The deviation from exponential recovery may well 
have been larger in the period just following the 
glitch, however, the large uncertainties on the pe- 
riodogram measurements of v 10^^ Hz) pre- 
vent any firm conclusion. However, since the cor- 
responding pulse TOAs cannot be unambiguously 
phase-connected, it is likely that large variations 
in V and v occurred during the 240-day period 
between the glitch epoch and when we regained 
phase-coherence. It is also possible that a second, 
smaller glitch {/S.v/v < 10~^) occurred during this 
period. The observed deviation from the expo- 
nential recovery decreases as the glitch recovers. 
Thus, in the closing months of 2008, the pulsar 
was rotating very regularly again, similar to its 
pre-glitch behavior. 

Figure H] shows post- glitch phase-coherent mea- 
surements of i> (a subset of the i> measurements 
shown in the middle panel of Figure [2]). The pre- 
glitch measurements are excluded here for clar- 
ity. The solid line is the derivative of the glitch 
model fitted to the measurements, clearly show- 
ing that there is significant deviation. The overall 
effect of the glitch recovery on ;> is clear, however, 
from MJD 54100-54300 the measurements de- 
viate from the exponential recovery by ^ 0.15%. 
The effect of this anomalous change in is not di- 
rectly evident in the measurements of v (which are 
dominated by the exponential recovery) but does 
help explain why the exponential glitch fit is not a 
satisfactory description of the data, and is clear in 
the residuals of the glitch fit in the bottom panel 
of Figure [21 
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Fig. 4. — Top panel: Post-burst /> measurements 
of PSR J1846— 0258. Black crosses are produced 
from short phase-coherent measurements of v and 
while filled red circles are from overlapping seg- 
ments of data, but are otherwise produced in the 
same manner. Open blue circles are from weighted 
least-squares fit to periodogram measurements of 
f. See Section [3] for details. Note the significant 
deviation from a simple exponential glitch recov- 
ery evident in the coherent z> measurements in the 
interval MJD 54100-54300. The solid line is the 
derivative of the exponential recovery model fitted 
to the frequency data. Bottom panel: Residuals 
of the v measurements and derivative of the expo- 
nential glitch recovery model. 



4. Bursts and pulsed flux 



The glitch during PSR J1846-0258's outburst 
was accompanied by a major pulsed fiux enhance- 
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ment (jGavriil et al.l 120081 ) . In order to quantify 
the radiative properties of the source during the 
gUtch recovery, we extracted its pulsed flux us- 
ing all available RXTE observations. First, we 
generated separate event lists for each PCU in 
FIT^ format using the standard FTOOL^. We 
then filtered our event lists such that we only pre- 
served photons in the 2 - 20keV band, from the 
first Xenon layer. The photon arrival times were 
then barycented using the source position and the 
JPL DE200 solar system ephemeris. We folded 
the filtered barycented photon arrival times using 
the ephemeris determined in our phase- coherent 
timing analysis using 16 phase bins. Using the 
folded profiles, we calculated the RMS pulsed flux 
in each PCU usi n g the Fourier method described 
by IWoods et al. (2004) keeping only the contri- 
bution from the 1st harmonic given the source's 
roughly sinusoidal profile. Not all the observations 
were pointed at PSR J1846— 0258, therefore we 
corrected for the reduced efficiency in each PCU 
due to the offset pointing using the collimator re- 
sponse of each PCU and the instrument attitude 
files. Finally, we averaged the pulsed flux in each 
PCU weighted by the fractional exposure of each 
PCU. We excluded the contribution from PCU 
because of the loss of its propane layer and the 
numerous detector breakdown events. Our pulsed 
flux time series is presented in the bottom panel of 
Figure [2] The event lists for each PCU created for 
the pulsed flux analysis were binned into 31.25 ms 
lightcurves and were searched for bu rsts using the 



burst search algorithm introduced in lGavriil et al 
(j2002( ). No additional bursts were found other 



than the 5 reported in ICavriil et al. ( 2008[ ). 

We checked for a correlation between torque 
and pulsed flux by plotting the pulsed flux against 
the spin frequency derivative in log-log space, 
shown in Figure [5] All measurements of u before 
and after the glitch and bursts are shown, includ- 
ing 3 measurements of z> obtained from weighted 
least-squares fits to v measurements obtained from 
periodograms. These 3 v measurements have un- 
certainties ^2 — 3 orders of magnitude larger than 
those for v measurements from coherent timing. 

The plot shows a possible correlation between 
flux and torque when both parameters are ex- 
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Fig. 5. — Observed relationship between pulsed 
flux and the spin-down rate of PSR J1846-0258. 
The measurements of v from phase-coherent tim- 
ing solutions before and after the outburst are 
shown as crosses. The three v measurements made 
from least squares fits to periodogram measure- 
ments of V are shown as filled circles. The best-fit 
power-law index to all data is ^ 8, however, the 
correlation is constrained by a single measurement 
of made when the timing noise was extremely 
large. Excluding this single point results in no 
significant correlation. See Section |4] for details. 



• ^http: //f its . gsf c .nasa. gov| 

^http: //heasarc .gsf c .nasa. gov/docs/sof t«are/f tools/| 
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treme. Fitting a power law to the data gives a 
power law index of ~8. However, this value is 
constrained by a single measurement of v when 
the flux, and timing noise, are both very large. 
Excluding this point, there is no significant corre- 
lation between flux and v. The observed variations 
in V during the later phase connected period, are 
at the ~ 3% level, while the variation in v imme- 
diately following the glitch are at the ~ 15% level 
(Figure ^ . Correspondingly small fluctuations in 
pulsed flux are not detectable in these data. 

5. Discussion 

5.1. Glitch properties 

Four X-ray bursts in PSR J1846-0258 coin- 
cided with the onset of a flux flare on 2006 May 
31 (MJD 53886). The pulsed flux decayed over 
~2 months and reached quiescence around the 
time of the fifth bur s t on 2006 July 27 (MJD 



53943) (jGavriil et a l] |2008|). Significant spec 



tral changes also occurred (iGavriil et al.l 120081: 



iKumar fc Safi-Harbl2008l : lKiiiDer fc Hermsenll2009h 
and flu x enhancement up to 30 keV was ob- 
served ( Kuiper fc Hermsenl 20091 ). Contempo- 
raneous with the sudden change in the X-ray 
emission from PSR Jf846— 0258, we observed a 
large glitch with an initial frequency increase 
oi V ^ 1.24(4f) X 10"^ The glitch decayed 
over f 27 (5) days, with a recovery fraction of 
Q — 8.7(2.5), resulting in a net decrease of the 
pulse frequency of Az^ — —9.52(9) x f 0~^ Hz. Fur- 
thermore, the timing behavior during the period of 
recovery is not well modeled by a simple exponen- 
tial function and measurements of v in particular 
are suggestive of a high level of timing noise for 
several hundred days following the glitch. 

This glitch and subsequent recovery reinforces 
that PSR Jf 846— 0258 underwent a period of mag- 
netic activity in 2006. This glitch is entirely 
different from the previous glitch in this source, 
which was radiatively silent, small in magnitude 
[lS.v/v = 2.5(2) X 10~^), dominated by a change 
in V (Ai>/i' = 9.3(f) x fO"**), and had no mea- 
surable recovery (LKGK06). This small glitch 
is si milar to those obse rved in the Crab pulsar 
(e.g., IWong et al.l l200l[ ) and other very young 
rotation-powered pulsars such as PSR B0540— 69 
(jLivingstone et al.l 120051 ) . It is unusual to have 
two such disparate initial /S.v magnitudes in the 



same source, particularly in such a young pul- 
sar, though this has been seen in some older pul- 
sars such as PSR B1737— 30, which has glitch 
magnitudes spanning four orders of magn itude 
( Lvne et al. 2000t Janssen &: Stappers 2006). In 
fact, the glitch reported here is the largest glitch 
ever observed in any of the pulsars with charac- 
teristic ages less than ^ 2 kyr (the Crab pulsar, 
B0540-69, B1509-59, and J1119-6127), none of 
which have experienced glitches with fractional 
magnitudes larger than /S.v/v ~ fO~*. 

The glitch reported here has a recovery frac- 
tion ofQ = 8.7±2.5. The measured value of 
Q > f implies that the net frequency change after 
the glitch recovery is negative, /S.Vp = —9.52(9) x 
10~^, as shown in Figure [2j A similar effect 
was recently observed in the AXP 4U 0142-1-61, 
though with much smaller ma gnitude of Q 



1.07 ± 0.02 (|Gavriil et al.ll2009l ). The negative 



change in v resulting from the over-recovery of the 
PSR J1846— 0258 glitch is similar in magnitude to 
the unresolved timing event seen in the magnetar 
SGR 1900-1- 14 in the month s before the giant flare 
in 1998 ( Woods et al. 19991 ). In that case, an en- 
hanced spin-down of the magnetar was observed 
over ^3 months. Well spaced timing observations 
around the time of the timing event were not avail- 
able, so no glitch could be resolved, if indeed one 
occurred. 

Thompson et al. ( 2000l ) attribute the observed 



behavior in SGR 1900-1-14 to one of two possi- 
bilities. The first is an increase in the magni- 
tude of r> by a factor of ^ 2.3, persisting for 
~ 80 days. The second is that a negative glitch, 
that is, a sudden spin-down occurred, with mag- 
nitude /S.V /v ~ 1 X 10^**. However, a timing event 
similar to that observed in PSR J1846— 0258 could 
also describe the data. It is curious, however, 
that such similar fractional changes in v should 
occur in two sources that experienced such dis- 
parate radiative changes, with the energy output 
from SGR 1900-1-14 being several orders of magni- 
tude larger than from PSR J1846-0258. 

AXP IE 2259-1-586 experienced a glitch con- 
temporaneous with 80 X-ray bursts, a fiux fiare 



and p ulse profile change s in 2002 (jKaspi et al 



20031 IWoods et al.l 12004^ . Neither this glitch 



nor the large PSR J1846-0258 glitch can be de- 
scribed with a single exponential decay model. 
The addition of an exponential growth compo- 
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nent better describes the IE 2259+586 glitch, 
however, the addition of a similar component for 
the PSR J1846-0258 glitch does not provide a 
significant improvement to the fit. Interestingly, 
the flux enhancement observed in IE 225 9+586 

20081) 



Zhu et al 



lasted much longer (> 2yr; 
than the glitch recovery time scale (r^ ^ 16 days) 
whereas the reverse is true for PSR J1846— 0258 
with Td = 127 days and a flux decay t imescale 



of 55.5 ± 5.7 days (|Gavriil et al.l l2008l) . The 



IE 2259+586 event can also be distinguished from 
the PSR J1846-0258 ghtch in that its recovery 
fraction is much smaller, with Q ~ 0.19. In 2001, 
the AXP IRXS J170849.0-400910 also experi- 
enced a glitch with recovery that was not well 
described by a simple exponential, and not im- 
proved wi th the addition of a s econd exponen- 
tial term Jkaspi fc Gavriill l2003t iDafl'Osso et al 



2003l : lDib et al.boOSi ). IWoods et all (|2004l ) argued 
that it is unlikely that IRXS J17Q849.0-400910 
experienced bursts or a pulsed flux flare associ- 
ated with this glitch because a flux flare would 
have had to decay on a time scale less than the 
glitch decay time scale. However, this is exactly 
the behavior observed from PSR J1846-0258, al- 
beit with a much longer time scale. Long-term 
spectral changes and flux variations have been 



claimed in IRXS J170849.0-400910 (Rea et al 



20051 : iCampana et al.ll2007t llsrael et al 



2007bh . 



An interesting characteristic of some AXP 
glitches is a period of enhanced spin-down im- 
mediately following the glitch, as observed in 
IE 1841-045, IRXS J170849 0-400910, and 
IE 2 259+586 (|Kaspi fc Gavriill l2003t IPib et al 



2008 ) . The additional instantaneous spin-down at 
the glitch epoch owing to the exponential recovery 
can be quantified as !>inst = —^i^d/Td (found by 
taking the derivative of the exponential term in 
Eq. 1 and setting t = tg). For the IE 2259+586 
glitch, i>i„st = (8.2±0.6)z >, while a typi cal value for 
a RPP is ~ 0.005i' (see IPeraltalliooel . and refer- 
ences contained therein). For PSR J1846-0258 
the instantaneous spin-down is v-mst — 0.15/>, 
larger than for any RPP glitch, but not as large 
as those measured for AXPs. It should be noted, 
however, that an enhanced spin-down is not ob- 
served in every AXP glitch. 

In PSR J1846-0258, there is evidence for a 
~200-day interval where v deviates significantly 
from an exponential glitch recovery, and it is possi- 



ble that further significant deviation occurred dur- 
ing the 240-day period of unconnected data di- 
rectly following the glitch. Perhaps v is varying 
in a stochastic fashion similar to that observed in 



Dib et al 



the AXP IE 10 48.1-5937 (jGavriil fc Kaspill2004 : 
2009t ). In this AXP, a ~year-long pe- 



riod of rapid v variations followed a large pulsed 
flux flare and a possible glitch in 2002. Another 
glitch in 2007 coincided with the onset of a pulsed- 
flux flare, again followed by stochastic variations 
in v. Another similarity between these sources is 
that the flux enhancement decayed away long be- 
fore the timing variations subsided. Alternatively, 
perhaps the variations are more simply attributed 
to timing noise, as is seen in many young RPPs. 
The behavior remains unusual however, since this 
is qualitatively very different from the mild timing 
noise observed in PSR J1846— 0258 prior to mag- 
netic activity, and such a large change in timing 
noise behavior is unprecedented among RPPs. 

A simple estimate of the transfer of rotational 
kinetic energy at the time of the glitch can be ob- 
tained by treating the star as a solid rotating body. 
Assuming the canonical neutron star moment of 
inertia of / = lO'^^gcni'^, the energy deposited in 
the solid crust at the time of the glitch owing to the 
increase in v is Ai? ~ (27r)^/i^Ai/ ^ 2 x lO'*^ ergs. 
However, it is well established that neutron stars 
are comprised of at least two components, a solid 
crust coupled to the core and a loosely coupled 
superfluid component in the crust. Taking the 
two-component model into account but not mak- 
ing any further assumptions about the nature of 
the two components, we can calculate the energy 
transferred between the components at the time 
of the glitch. The first constraint is that angular 
momentum is conserved, that is, 



Ic^v = IsF^l'SF, 



(2) 



where Ai/ is the observed increase in spin- 
frequency, and l^usF is the unknown change in 
spin-frequency of the superfluid, and Ic and Isf 
are the moment of inertia of the solid crust and 
core, and superfluid, respectively. An estimate 
of the energy contained in the glitch can then be 
calculated as 

Eg = AEc - AEsF - (2^)2/cAj.(i/ - i^sf)- (3) 

The frequency lag between the crust and super- 
fluid can be estimated from the glitch decay time. 
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Td, as 



Tr 



(4) 



where Tc is the characteri stic age of the pulsar 
(|Shapiro fc Teukolskvll983h . For PSR J1846-0258, 
the frequency lag \s vsf — ^ 0.0013Hz, giving 
a glitch energy of 7 x 10'^^ ergs for the measured 
glitch parameters (see Table [Ij . If the measured 
glitch decay timescale, Td, in this case is domi- 
nated by processes external to the neutron star (as 
discussed in Section [O)) . any other glitch decay 
timescale would occur on a shorter timescale than 
the one observed, so this estimate of the frequency 
lag and thus glitch energy can be considered an 
upper limit. 



For a distance to the pulsar of 6 kpc ([Leahy fc Tian 



20081 ) the energy estimated to have been released 
in the bursts and flux flare, assuming isotropic 
emission, is (3.8 - 4.8) x 10'*i(d/6kpc)2 ergs (2 - 
60keV). A new est imate of the distance to the 
pulsar of '^10 kpc ( Su et al.l 12009) increases the 
amount of energy contained in the radiative out- 
burst to (1.1 - 1.3) X 10*2 (d/10kpc)2 ergs (2 - 
60keV). For either distance, the energy contained 
in the radiative outburst is several orders of mag- 
nitude larger than that contained in the glitch. 
This is similar to the bursts and glitch from 
IE 2259-1-586, for which the energy contained in 
the glitch was ^2 orders of magnitude less than 
the energy containe d in the bursts and flux flare 



( Woods et al. l2004l) . This is suggestive that the 



glitch alone is not responsible for the radiative 
ou tburst, in contrast to the a rgument put forward 
bv lKuiper fc Hermsen ( 20091 ). Moreover, there is 
no evidence for this event (or in any other similar 
AXP event) that the glitch preceded the radiative 
event. This could be tested only by sensitive con- 
tinuous X-ray monitoring of this and other similar 
sources. 

5.2. Physical models for "magnetic glitches" 

Rotation-powered pulsar glitches are thought 
to arise from differential rotation in the neutron 
star, where the crust contains superfluid neutrons 
rotating more rapidly than the surroun ding mat- 



ter (e.g., lAlpar et al.lll984l : lAlpar fc Pin es 1993) 



The angular momentum of a rotating superfluid 
is quantized in vortices which are thought to be- 
come pinned to nuclei in the star's crust. The 
vortex lines are therefore under extreme stresses 



due to the differential rotation between the crust 
and the superfluid. It is thought that the vor- 
tex lines experience sudden unpinning, resulting 
in the transfer of angular momentum to the crust, 
observed as an increase in the spin-frequency of 
the neutron star. Magnetar glitches may instead 
be triggered by strong internal magnetic fields as 
the crust is deformed, either plasticall y or cracked 
violently ( Thompson fc DuncanI [l996l ) . This idea 



is supported by the large number of glitches now 
observed to occur at the same epoch as magneti- 



flares (e.g., 


Kasoi et al. 


Ilsrael et al. 


l2007ah. 



20031 iDib et al.l l2009t 



The physics underlying glitches with Q > 1 is 
unclear. The classical glitch model of vortex un- 
pinning in the superfluid crust of the neutron star 
does not readily produce such dramatic glitch re- 
coveries. One possibility is that some parts of the 
superfluid are in fact rotating more slowly than 
the crust. Then the initial v increase would be 
from a transfer of angular momentum from a more 
rapidly rotating region of the superfluid to the 
crust, which is then followed by a transfer of angu- 
lar momentum from the crust to the more sluggish 
region of the s uperfluid. This i s the explanation 
put forward by [Thompson et al.l ((2000, ) to explain 
the net spin-down event in SGR 1900-1-14. They 
argue that regions of slowly rotating superfluid can 
occur in magnetars because vortex motion is domi- 
nated by advection across the neutron star surface 
by the deforming crust and that gradual plastic de- 
formation of the neutron star crust will cause the 
superfluid to rotate more slowly than the crust. 
However, it is not clear whether such behavior is 
expected in a neutron star with a magnetic field of 
~ 5 X 10^'^ G, spinning relatively rapidly compared 
to the magnetars. 

That the recovery of the PSR J1846-0258 
glitch so far overshoots the initial frequency in- 
crease is suggestive of an external torque follow- 
ing the glitch. Previous evidence for an external 
torque includes the large fraction of / implied to 
have decoupled in the IE 2259-1-586 glitch. How- 
ever, this may also be explained with a core glitch, 
that is, where the core of the pulsar (and thus a 
large f raction of I) decouples temporarily from th e 
crust (|Kaspi fc Gavriilll2Q03l: IWoods et al.ll2004[) . 
In PSR J1846-0258, however, the post-glitch re- 
laxation amplitude is much greater than the initial 
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glitch amplitude, offering support for the idea that 
the post-glitch spin-down behavior results from an 
external source. 

One possibility is that a magnetic field twist 
responsible for the X-ray bursts and flux en- 
hancement also affects the spin-down of the pul- 
sar. In this case the observed recovery is driven 
by the propagation of magnetic field untwisting 
(similar to a shock wave) through the magne- 
tosphere. During this process, the spin-down 
of the star may increase because the effective 
magnetic field has increased. When the "shock" 
reaches the light cylinder, which can take place 
on few month time scales, t he spin-do wn should 
return to its pre-burst value (jBeloborodoy 2009). 
This theory also allows for non-monotonic behav- 
ior in the spin-down after an ev ent, as observed 
in both AXP IE 10 48.1-5937 (|Gavriil fc Kaspi 
20041: iDib et "aDl2Q09l) and PSR J1846-0258. This 
model allows for a delay between flux variations 
and the onset of timing variability, as was ob- 
served in IE 1048.1—5937, though not observed 
in PSR J1846-0258. However, as in the model of 
variably rotating superfiuid, it is not immediately 
clear how the more rapid rotation and smaller 
magnetic field affect the relevance of this model. 
In particular, the different light cylinder radii of 
PSR J1846-0258 compared to the AXPs should 
result in different shock propagation times. 

Alternatively, it has been proposed that fall- 
back disks from the supernova explosion creat- 
ing the neutron star could be interacting with 
magnetars and be respo nsible for some of the 



Chatteriee et all l2000t 



observed em ission (e.g. 
Alpaiil200l[ ). In this case, the initial X-ray bursts 
could irradiate a fossil disk commencing a period 
of disk activity. The interaction between the neu- 
tron star and a disk could cause the enhanced 
spin-down, which would decay as the disk cooled. 
However, in the framework of this model, it is dif- 
ficult to understand how accretion causing varia- 
tions in the spin-down rate could continue for so 
much longer than the pulsed- fiux enhancement. 

5.3. Magnetar and High-B Radio Pulsar 
Properties 

Another RPP, PSR J1119-6127, has very sim- 
ilar spin properties to PSR J1846-0258. Its P, E, 
and Tc are all similar, and notably, it has a sim- 
ilarly large magnetic field of i? = 4.1 x 10^^ G 



(jCamilo et al.l I2OOOI) . This pulsar has shown 
some indication of unusual X-ray emission. No 
magnetospheric X-ray emission is detected, but 
thermal pulsations with a ^ 75% pulsed frac- 
tion and a la rge surface temperature are detected 
(jGonzalez et al. 2005). No direct evidence of mag- 
netic activity (i.e. bursts or fiux enhancements) 
is present in this pulsar. Given the similarities 
between PSR J1119-6127 and PSR J1846-0258, 
both sources, as well as other high B-field RPPs 
are currently being monitored with RXTE for sim- 
ilar magnetic activity. 

PSR J1846-0258 may be related to the tran- 
sient AXPs (TAXPs). PSR J1846-0258 appears 
to be a typical RPP for > 95% of the time, 
with brief periods of magnetic activity occur- 
ring approximately once a decade. The TAXP 
XTE J1810— 197 increased in brightness by a fac- 
tor of ^100 and was subsequently visible for sev- 



eral years as a pulsed X-ray source (llbrahim et al 



_2004; Halpern_& Gotthch 2005). Another TAXP, 
IE 1547—5408 was detected as a pulsed radio 
source after an X-ray outburst in which the fiux 



incre ased by at least a factor of 16 (jCamilo et al 
2007t ). Interestingly, both XTE J1810-197 and 
IE 1547-5408, two bona fide TAXPs, are the 
only two magnetars with detected radio pulsa- 
tions. By contrast, no radio pulsations have been 
detected from PSR J1846— 0258 despite extensive 
sear ches both before and after the magnetic activ- 
ity (|Kaspi et al.lll996HArchibald et al.il2008[) . 



6. Conclusions 

In RXTE observations of PSR J1846-0258, we 
have observed a large glitch with an unusual quasi- 
exponential over-recovery of v and substantial 
timing noise contemporaneous with X-ray bursts 
and a fiux increase. These observations strengthen 
the tie between magnetic activity in neutron stars 
and unusual glitc h activity, as h as been previ- 



ously noted (e.g., iDib et al.ll2008l ). A glitch with 



recovery fraction Q > 1 has never before been 
observed from a rotation-powered pulsar and is 
not compatible with the standard model of pul- 
sar glitches. The unusually large glitch recov- 
ery reported here for PSR J1846-0258, as weU 
as the radiative changes occurring contemporane- 
ously with the PSR J1846-0258 glitch and several 
AXP glitches, together provide the best evidence 
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that there are physical differences between typical 
RPP glitches and some glitches observed in mag- 
netars. 

Ongoing timing observations of PSR J1846— 0258 
are required to obtain a deterministic post- 
outburst braking index measurement that is un- 
contaminated by long-term ghtch recovery. A 
measurement of the braking index for this source 
after outburst is of considerable interest as a 
change would strongly suggest that a magnetic 
reconfiguration occurred at the time of the out- 
burst. 
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Table 1 



Glitch Parameters for PSR J1846-0258. 
Parameter Value 


Epoch (Modified Julian Day) 


53883.0(3.0) 




4.0(1.3) X 10"^ 




0.0041(2) 


Td (days) 


127(5) 


Q 


8.7(2.5) 




-9.52(9) X 10-^ 
10.8(4) X 10-5 



Figures in parentheses are uncertainties in the last digits quoted and are the estimated Itr uncertainties 
from the fitted exponential glitch recovery model using Ax^ = 1 contours. The model is fitted to data with 
uncertainties that are increased by a factor until ~ 1 (see Section [3] for details). 
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